MARKOV CHAIN APPROXIMATIONS FOR TRANSITION DENSITIES OF 

LEVY PROCESSES 



ALEKSANDAR MIJATOVIC, MATIJA VIDMAR, AND SAUL JACKA 

Abstract. We consider the convergence of a continuous-time Markov chain approximation X h , 
h > 0, to an K-valued Levy process X. The state space of X h is an equidistant lattice and its 
Q-matrix is chosen to approximate the generator of X. Under a general sufficient condition for 
the existence of transition densities of X, we establish sharp convergence rates of the normalised 
probability mass function of X h to the probability density function of X . 

1. Introduction 

Discretization schemes for stochastic processes are relevant both theoretically, as they shed light 
on the nature of the underlying stochasticity, and practically since they lend themselves well to 
numerical methods. Indeed, there has been a plethora of publications devoted to the subject, see 
e.g. [9], and with regard to the pricing of financial derivatives [7] and the references therein. 

Levy processes, in particular, constitute a rich and fundamental class with applications in diverse 
areas such as mathematical finance, risk management, insurance, queuing, storage and population 
genetics (see e.g. [10J). There is a wealth of existing literature concerning approximations of Levy 
processes in one form or another and a brief overview of simulation techniques is given by |13| . In 
continuous time |8j approximates by replacing the small jumps part by a diffusion, and discusses also 
rates of convergence for Eji/oX^], where g is real- valued and satisfies certain integrability conditions, 
T is a fixed time and X the process under approximation; [5] approximates by a combination of 
Brownian motion and sums of compound Poisson processes with two-sided exponential densities. In 
discrete time, Markov chains have been used to approximate the much larger class of Feller processes 
and [3] proves convergence in law of such an approximation in the Skorohod space of cadlag paths, 
but does not discuss rates of convergence; [T3] has a finite state-space path approximation and 
applies this to option pricing together with a discussion of the rates of convergence for the prices. 

This paper studies the rate of convergence of a weak approximation of a Levy process X by 
a continuous-time Markov chain (CTMC). Our main result, Theorem [2j establishes the precise 
convergence rates of the normalised probability mass functions of the approximating Markov chains 
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to the transition densities of the Levy process for the two proposed discretisation schemes, one in 
the case where X has a non-trivial diffusion component and one when it does not. More precisely, 
in both cases we approximate X by a CTMC X h with state space := KL and Q-matrix defined 
as a natural discretised version of the generator of X. This makes the CTMC X h into a continuous- 
time random walk, which is skip-free (i.e. simple) if X is without jumps (i.e. Brownian motion 
with drift). The quantity 

n(h) := / |x|dA(a;), 
J[-l,l]\[-h,h] 

where A is the Levy measure of X, is related to the activity of the small jumps of X and plays 
a crucial role in the rate of convergence. We assume that either the diffusion component of X is 
present (a 2 > 0) or the jump activity of X is sufficient (Orey's condition [14, Prop. 28.3], see also 
Assumption [T] below) to ensure that X admits continuous transition densities pt : T(x,y) (from x 



at time t to y at time T > t). Let P£ T (x,y) := P(X£ = y\X[ 
transition probabilities of X h and let 

A T -t(h) 



denote the corresponding 



sup 

x,yeZ h 



Pt,r(x,y) - -P ttT (x,y) 



The following table summarizes our result (here for functions / > and g > 0, we write / = 0(g) 
for limsup/^o f(h)/g(h) < oo — if g converges to 0, we say / decays no slower than g): 



A T -t(h) 


a 2 >0 


a 2 = 


A(R) = 


0(h 2 ) 


X 


< A(R) < oo 


O(h) 


X 


A(M) = oo 


0(hn(h/2)) 



We also prove that the rates stated here are sharp in the sense that there exist Levy processes for 
which convergence is no better than stated. 

The rate of convergence depends on the Levy measure A only in the case when X has infinite 
activity jumps. If X is a compound Poisson process (i.e. A(R) G (0, oo)) without a diffusion 
component, there is always an atom at zero and hence the convergence is studied only in the case 
a 2 > 0. In the infinite activity case, the definition of k implies k(Ji) = o(l/h) (for functions / > 
and g > 0, we use the notation / = o(g) to mean lim/40 f(h)/g(h) = — if converges to 0, we 
say / decays strictly faster than g). However the convergence of hK,(h/2) to zero, as h 1 0, can be 
arbitrarily slow. 

We point out that a result, which is less explicit, but similar in spirit to ours, is to be found in 
PQ for solutions to SDEs: for the case of the Euler approximation scheme, the authors there also 
study the rate of convergence of the transition densities. 

The proof of Theorem [2] is in two steps: we first establish the convergence rates of the charac- 



teristic exponents (Subsection 3.2) of X t and X^. In the second step we apply these to study the 
convergence of the transition densities (Section [4| via their spectral representations (established in 



Subsection 3.1). Note that in general the rates of convergence of the characteristic functions do 
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not carry over directly to the distribution functions. We are able to follow through the above pro- 
gramme by exploiting the special structure of the infinitely divisible distributions in what amounts 
to a detailed comparison of the transition kernels Pt,T(%,y) and P^ T {x,y). 

By way of example, note that if A([— 1, 1]\[— h, h\) ~ l/h l+a for some a G (0, 1), then n(h) ~ h~ a 
and the convergence of the normalized probability mass function to the transition density is by 
Theorem [2] of order since k(0) = 00 and Orey's condition is satisfied (we use the notation 

/ ~ g to mean lim^o f(h)/g(h) G (0, 00)). In particular, in the case of the CGMY [I] (tempered 
stable) or /3-stable [TJl p. 80] processes with stability parameter j3 G (1, 2), we have a = (3 — 1 and 
hence the convergence of order h 2 ^ 13 . 

Observe also that in terms of pricing theory, the probability density function of a process can be 
viewed as the Arrow-Debreu state price, i.e. the current value of an option whose payoff equals the 
Dirac delta function. The singular nature of this payoff makes it hard, particularly in the presence 
of jumps, to study the convergence of the prices under the discretised process to its continuous 
counterpart. Theorem [2] can be viewed as generalisation of such convergence results for the well- 
known discretisation of the Black-Scholes model (see e.g. [11J for the case of Brownian motion with 
drift). 

The rest of the paper is organised as follows. Section [2] introduces the setting by specifying the 
Markov generator of X h , and precisely states the main result. In section [3] we find integral expres- 
sions for the transition kernels by applying spectral theory to the generator of the approximating 
chain. In section [4] this allows us to establish convergence rates for the transition densities. 

2. Definitions, notation and statement of results 

2.1. Setting. Let *(p) := -\a 2 p 2 + ifip + J R (e ipx - 1 - ipxl[_ V}V] (x))d\(x) (p G M) be the 
characteristic exponent corresponding to a Levy process X with characteristic triplet (a 2 , A,/i)g — 
c(y) := l[_y;y](y) is the cut-off function and V is 1 or 0, the latter only if J^_ 1 ^ \x\d\(x) < 00. 
Then X is a Markov process with transition function Pt j<(x, B) := P(Xx-t G B — x) (0 < t < T, 
x G R and B G B(R)) and <f>x t ip) = E[e ipXt ] = exp{i^(p)}. We refer to [2], H3] for the background 
on Levy processes. 

Fix h > 0. Consider a CTMC X h = {X^)t>o approximating our Levy process X (in law). We 
describe X h as having [12] state space := KL := {hk : k G Z}, initial state Xq = 0, a.s. and an 
infinitesimal generator L h given by a spatially homogeneous Q-matrix Q h (i.e. , depends only 
on s — s', for {s, s'} C Z/J. Thus L h is a mapping defined on the set Z°°(Z^) of bounded functions 
/ on Z h , and L h f(s) = Zs'eS Q^fW)- 

It remains to specify Q h . To this end we discretise on the infinitesimal generator L of the 
Levy process X, thus obtaining L h . Recall that [IU p. 208]: 

Lf(x) = a -f"(x)+ f ,f'(x)+ / [f{ X + y)-f{ X )-yf'( X )l [ _ vv] (y)]d\{y) 
1 JR 



MARKOV CHAIN APPROXIMATIONS FOR TRANSITION DENSITIES OF LEVY PROCESSES 



1 



for i6K and every twice continuously differentiable / : M — > K with /, /', /" vanishing at infinity. 
We thus have discretisation scheme 1, for s G 7Lh and / : — >■ K vanishing at infinity: 

-K t< „, 1 / 2 , K\ f(s + h)+f( S -h)-2f(s) , / ,X /( S + / l )-/( S -/ l ) 



a) + £ [/(-+ ^ -/w] 4 

s'ez h \{o} 

where we have introduced := A(^4g) with := (s — h/2,s + h/2] if s > and := [s — h/2,s + 
h/2) if s < (s G Z/,\{0}); finally (letting := [-h/2,h/2]; 5° G {0, 1}): 

c h :=5° [ y 2 l { - Vy V]{y)d\{y) £ s l h-V^WKv)- 

I A h I Ah 

JA o s&Z h \{0} JA » 

We reserve the right to set Cq to even if V 7^ 0, hence the presence of 5° G {0, 1}. 
Note that Q h has nonnegative off-diagonal entries for all h for which: 

and in that case Q h is a genuine Q-matrix (due to spatial homogeneity, its entries are then also 
uniformly bounded in absolute value). When a 2 > 0, it will be shown that Q always obtains, at 



least for all sufficiently small h (see Proposition 11 ), but not in general. When a = we shall thus 



employ an alternative discretisation for which the above condition holds vacuously. The changes 
are as follows. 

First, (|T| now reads (letting sgn := l(o,oo) ~~ l(-oo,o)) : 

lVW = \{^ + 4) /(' + »> + /fc - M -VM + E [/(s + s . + ftsgn(s0) _ /(s)]4 + 



(3) + [p-fi ) 7 l[0,oo)(A* - M J + (M - M J 7 1 (-oo,0] (A 1 _ M ) 



s'GZ h \{0} 

>)z 



So, in place of (/x — fi h ) ^ s+h ^~^ s ~ h ^> we have rather (fi — /U ' l )Zi£±^zZM or (^ _ ^h^ f(s)-f(s-h) 
accordingly as [i — [i h > or \i — fi h < 0. We have also modified slightly the final summation in L h . 
This will turn out to be advantageous when considering the convergence of the transition kernel 



with a = (see the estimate in (25)) and is made solely for this purpose: it is possible that one 
could establish convergence rates even without this modification but our proof technique relies on 
it. 

Second, we have: 



s&Z h \{0} 



(s + fasgn(fl)) / h-vy\(v)d\(y). 



Everything else remains unchanged. We shall refer to this as discretisation scheme 2. Most of 
our arguments will be the same under both and when it is not so, we shall make this explicit. 
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2.2. Summary of results. The assumption alluded to in the introduction is the following (we 
state it explicitly when it is being used): 



Assumption 1. Either a 2 > or Orey's condition holds: 



3eG (0,2) 



such that 



lim inf 



1 



,2-e 



u 2 d\{u) > 0. 



riO r* - J[-r t r] 

The usage of the two schemes and the specification of 5° and V is as summarized in Table [l] - 
by contrast to Assumption [T] we maintain this table as being in effect throughout this subsection. 

Table 1 . Usage of the two schemes and of 5° and V depending on the nature of a 2 
and A. 



Levy measure/diffusion part 


a 2 > (scheme 1) 


a 2 = (scheme 2) 


A(R) < oo 


V = 6° = 


X 


k(0) < oo = A(M) 


V = 1 and 5° = 


k(0) = oo 


V = 5° = 1 



Under Assumption [T] for every t > 0, <j>x t £ L 1 (m) where m is Lebesgue measure and (for 
< t < T, {x,y} C M): 



(4) 



Pt,r(x,y) = — I exp{ip(x - y)}exp{*(p)(T -t)}dp. 



Similarly, with ty h denoting the characteristic exponent of the compound Poisson process X h (for 
< t < T, y G Zft, P X h-a.s. in x E Z^): 



(5) 



ip/; T (x, y) = Y j\ exp{ip(x - y)} exp{V h (p)(T - t)}dp. 



Note that the RHS is defined even if P(X£ = x) 
so. Finally, for < t < T and x,y £ Z^ define: 



and we let the LHS be this value when this is 



(6) A r -*(/i):= sup D$ T (x,y) where D% T (x,y) 
{x,y}cz h 

The main result of the paper can now be stated. 



Pt,r(x,y) - ^Pt )T {x,y) 



Theorem 2 (Convergence of transition kernels). Under Assumption [i| whenever s > 0, the con- 
vergence of A s (h) is as summarized in the following table. In general convergence is no better than 
stipulated. 





A(M) = 


< A(M) < oo 


k(0) < oo = A(R) 


k{0) = oo 


a 2 >0 


0(h 2 ) 


O(h) 


0(h) 


0(hn(h/2) 


a 2 = 


X 


X 
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Remark 3 (Convergence in distribution). X h converges to X weakly in finite- dimensional distri- 
butions as h ], 0. 

3. Transition kernels and convergence of the characteristic exponents 

3.1. Integral representations. First we note the following result (its proof is essentially by the 
standard inversion theorem, see also [14, p. 190]). 

Remark 4. Under Assumption [IJ Xt (t > 0) admits the continuous density fx t {v) = 
5tF Jr e ~ ipV( l ) Xt{p)dp- For some {P, C, e} C (0, oo) depending only on {A, a 2 } and allp G M\(—P, P): 
\4>x t (p)\ — exj p{~Ct\p\ 6 } . When a 2 > 0, e = 2 and otherwise e takes the value in Orey's condition 
in Assumption [7J In particular, the law P^t{x, •) is given by Q). 

Second, to obtain ^ we apply some classical theory of Hilbert spaces, see e.g. [6]. 

-£,£] -> C be given by g ' - 1 u 
constitute an orthonormal basis of the Hilbert space L 2 ([— ? , ?]). 



Definition 5. For s £ Z h let g s : [-f , f ] -» C be given by g s (p) := \j ^e lsp . The (g s ) se z h 

'A> hi 



Let A G Z 2 (Z/J. We define: ThA := EseZ ^( s )^ s- ^ e ^nuerse of this transform T h x : 
L2 (hM)^ 2 ( Z fc) is 9iven by: 

= {(/>,g s ) ■= / ^ftrfm 

for<f>€ L 2 {[-l,l\) andseZ h . 

Definition 6. For a bounded linear operator A : Z 2 (Z^) — > Z 2 (Z^) ; we say Fa ■ [—ir/h,Tr/h] — > R is 
its diagonalization, if T^ATf^^ = Fa^ for all (f) G L 2 ([— ^]). 



We now diagonalize I/' 4 , which allows us to establish Q. The straightforward proof is left to the 
reader. 



Proposition 7. Fix C G I (Z/j). T/ie following introduces a number of bounded linear operators 
A : l 2 (Zh) — )• Z 2 (Z/J and giwes i/ieir diagonalization. With f G l 2 (Zh), s G Z^, p G [— ? , ?]: 
^ A,/( S ) := F A » = 2^Jfti. 

(H) V ft /(a) := /(s+fe) ~ /(s ~ fe) . F v » = i?^. f/nder sc/ieme jg we Ze* V+/(s) := 

fre S p. Vfc/(a) := /(s) -{ (s - fe) ; and then F^p) = fresp. F v -(p) = ±=^J. 

(m) Lcf(s) := E a 'gz h (/(« + " /00)CV)- Fl c (p) = E seZfc C( S )(e^ - 1). f/nder 
scheme 2 we let Lcf(s) := Es'ez (f( s + s ' + tegn(s')) — f(s))C(s') and then Fl c (p) := 
E se z^(«)(e i{s+/lsgn(s))p -l)- 

As A is finite outside any neighborhood of 0, l/ 1 ^^) (as in Q, resp. ([3])) is a bounded linear 
mapping. We denote this restriction by L h also. Its diagonalization is then given by ^ h := F L h, 
where, under scheme 1, 
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sez h \{o} 



and under scheme 2, 



i/ip _ -i 1 _ -i/ip 
VHP) = (/X-/)^^l [0 ,oo)^-/) + (^-/) ^ l(-cx>,0](M-/) 

(8) + (a 2 + cg) (C0S(/ ^ 2 ) ~ 1) + ^ c ^(>+^gn(s)) P - l) 

sez h \{0} 

(with p 6 [— but we can and will view ty h as defined for all real p by the formulae above). 
Under either scheme, is bounded and continuous as the final sum converges absolutely and 
uniformly. 

Lemma 8. For scheme 1 under |Sp and always for scheme 2, for every < t < T, y £ Z/j and 
Pjf/j -a.s. in x £ 7;h |5p obtains, i.e.: 

— 

P(x£ = y\X? = x) = — f h e W {ip(x-y)}exp{^ h (p)(T-t)}dp. 

Proof. Note that: P(X£ = w|X t h = s) = (e( T -*) Lh !{,,}) (a;). Thus ^ follows directly from the 
relation ThL h TZ = ^> h - (where \E fft- - is the operator that multiplies functions point-wise 

In what follows we study the convergence of @ to (|]) as h I 0. These expressions are particularly 
suited to such an analysis, not least of all because the spatial and temporal components are nicely 
separated and homogeneous. 

One also checks that for every t > 0: 

</> X H(p) = E[J* x ?]=exp{m h ( P )} 
(hence X h are compound Poisson processes [144 P- 18, Definition 4.2]). 

3.2. Convergence of the characteristic exponents. We introduce for 

cos(hp) - 1 p 2 
Up) ■= p + y 

and, under scheme 1: 

. / sm(hp) 
9h{P> ■= i I — ^ P 

l h (p) := c ^os(M-l _^sinM + ^ c A (e ^_ 1 )_/'( e *«_i_ i p U i [ _ v;v](tt) ) (iA(l 

«6Z h \{0} • 7R 



MARKOV CHAIN APPROXIMATIONS FOR TRANSITION DENSITIES OF LEVY PROCESSES 



8 



respectively, under scheme 2: 



e 



ihp j j g—ihp 



ghip) ■= — r-^i[o,oo)0 - m' 1 ) h 7 — 1(-oo,o](m-m' 1 ) 



l[0,oo)0 - M ) + 7 l(-oo,0](M - P ! 



+ <£ ( e <(-+^W'»)P - x) - I (e ipu - 1 - iput^w^u)) dX(u). 



sez h \{o} 



Thus: 



y h - * = : + m + fo. 



Next, two elementary but key lemmas. The first (its proof, which relies on the Taylor expansion 
at and the Mean Value Theorem, is left to the reader) concerns some elementary trigonometric 
inequalities as well as the Lipschitz difference for the remainder of the exponential series fi(x) := 
Sfcli+i ^fr~ ( x G ^ S N U {0}). These estimates will be used again and again in what follows. 
The second lemma establishes key convergence properties relating to A and is proved below. 

2 4 3 

Lemma 9. For all real x: < cos(x) — 1 + ^- < ^j, < sgn(rr)(x — sin(x)) < sgn(a;)^ and 
< x 2 + 2(1 — cos(x)) — 2xsin(x) < x /4. Whenever {x,y} Cl« have (with 5 := y — x): 

(1) \e ix - e iy \ 2 = \e ix - 1 - (e iy - 1)| 2 < 5 2 . 

(2) \e lx - 1 - ix - (e iy - 1 - iy)\ 2 < 5 4 /4 + 8 2 x 2 + |<5| 3 |x| . 

(S) \e lx -l-ix+x 2 /2-(e i y-l-iy+y 2 /2)\ 2 < 5 6 /M+\5\ 5 \x\/6+(5/12)5 4 x 2 + \5\ 3 \x\ 3 /2+5 2 x* /4. 
Moreover, for all I G NU{0} ; one has \fi(x) - fi(y)\ < {V2) l2 \6\(\x\ l + \5\ l ). 

Lemma 10. For any Levy measure X one has for the two functions (given for 1 > 5 > 0): 
M (5) := 5 2 Jrj 1] \ ( _ (55) d\(x) and M x {5) := 5 J, 1 q\t SS s \x\dX{x) that M (5) -)• and Mx{5) -)• 
as 5 10. If, moreover, Jrj ^ |x|<iA(x) < oo, then S f,^ i]\rg$\ dX(x) — ^ as 6 10. 

Proof. Indeed let \x be the finite measure on ([—1, l],Br_i ji) given by fJ-(A) := J A x 2 dX(x) (A Borel 
subset of [-1, 1]) and let f$(x) := (f ) 2 1[_i i i]\(_ 5i(5 )(x) and f}(x) := |fy 1 [_ i , i] \(_^, 5) (^) be functions 
on [—1, 1]. Clearly < fg, fg<l and f$, fg—tO pointwise as 5 I 0. Hence by Lebesgue dominated 
convergence theorem (DCT), we have Mq (5) = f fgd/j, and M±(S) = J f l {5)dii converging to 
J Odfi = as 5 4- 0. The "finite first moment" case is similar. □ 



Proposition 11. Under scheme 1, with a 2 > 0, holds for all sufficiently small h. 
Proof. If V = this is immediate. If V = 1, then (via a triangle inequality) 

h\fx h \ < h 



s&Z h \{0} 



V sf ![_!,!] (y)dA(y) <h V / \s - u + u\t { _ ltl] {y)dX{y) 

■t, \ mi JAi v fnl J A 1 } 



s&Z h \{0} ' 



< h [ Ja([-1, 1]\[-/i/2, h/2]) + / |«|dA(«) ) ->• 

-l,l]\[/i/2,h/2] 
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as h 1 by Lemma 10 Eventually the expression is smaller than a 2 > and the claim follows. □ 

Furthermore, we have the following inequalities, which together imply an estimate for \^f h — 
^f\. For ease of notation in this proposition and the proofs to come, we introduce the notation 
(5 G (0,1]): ((5) := 5 J { _ im _ s>s] \x\dX(x), 7 (<5) := 8 2 J^^^ dX(x), c := A(M), b := «(0), 
d := A(M\[— 1, 1]). Recall also the definition of the sets Ag following 0. 

Proposition 12 (Convergence of characteristic exponents). For all p £ R: < fh(p) < p 4 h 2 /4! 
and < isgn(p)gfi(p) < /i 2 |p| 3 /3! (resp., under scheme 2, \gh(p)\ < hp 2 /2\). Moreover: 

(i) when c < 00; with 5° = V = 0: |//i(p)| < c\p\h/2 (resp. under scheme 2, \lh{p)\ 5; c M^j- 

(ii) whenb < 00 = c; wii/i 5° = and V = 1; for all h < 2: |^(p)| < | (|p|d + p 2 &) + (p 2 + |p| 3 )o(/i) 
(resp. under scheme 2, \lh{p)\ < §£> 2 &) + (p 2 + |p| 3 +P 4 )o(/i)J where o(h) depends 
only on A. 

(in,) whenb = 00; with 5° = V = I; for all h< 2: |Z ft (p)| < p 2 (({h/2) + ^(/i/ 2 )) +(bl +P 2 + M 3 + 
p 4 )0(h) (resp. under scheme 2, \l h (p)\ < p 2 [4,((h/2) + ^j(h/2)] + (|p| +p 2 + \p\ 3 + p 4 )0{h) ) 
and 0(h) depends only on A. Note here that we always have 7 < C an d that £ decays strictly 
slower than h, as h 1 0. 

Proof. The first two assertions are transparent by Lemma [9] — with the exception of the estimate 
under scheme 2, where (with 5 := hp): 

16 2 



\9h(p)\ 



2<fsin(<J) + 2(1 - cos(<5)) < 



hp z /2\. 



h v w v w/ ~ h 2 

For the final set of inequalities, the transposition of the estimates under scheme 1 to the ones under 
scheme 2 is essentially straightforward. Notable differences are given paranthetically in square 
brackets. 

If c < 00 (under 5° = V = 0): 



V c h s {e' sp -l)~ f (e' pu - l)dA(u) 

M[-ft,ft] 



s£Z h \{0} 



\ ^ h isp 

3£Z h \{0} 



" 2 ' 2 1 



E 

iS2 h \{0} 

E 



1-e 



ip(u — s) 



dA(w) 



d\(u) < \p\h\ K\ 



ft ft 
2' 2 



/2 



where in the second inequality we apply Lemma[9]and the first follows from the triangle inequalities. 
Finally | f<_ h/2>h/2] (e ipu - l)dX(u)\ < X([-h/2,h/2])\p\h/2 and the claim obtains. The reader will 
have no trouble reassuring himself that, under scheme 2, one obtains the bound as given. 
If b < 00 = c (under 5 = and V = 1; for all h < 2), denoting rj(S) := frg/n \x\dX(x): 



/( e ,pu — 1 — zptil [_i,i] (it) I d\(u) 
-ft, ft] ^ ' ' 



(9) 



(cos(pu) — 1) dX(u) 

■ ■ 41 

< ^(ft 2) • ( 



(sin(pu) — pu) dX(u) 



p\'h* V (h/2). 
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|p 9h(p)\ = 



h ( sin(/ip) 



< /i 2 |/||p| 3 /3! 



(10) 



i\dX(i 



!GZ h \{0} 

£ 

«ez fe \{o} 



2_] c s (e lsp — 1) — ipfJ. h — / (e lpu — 1 — ipul [_i i:L ] (u))dA(w) 

> x rni •/R\[-4.£l 



A h 



ipul [-1,1] («) + ipslf-i,!] (u) dA(it) 



dA(u) + p 



[-i,i]\[-f ,f] 



jujdA(u) + 



5 fh 



8 V 2 



dA(it) 



(11) 



, fc, , , . h 2, . 5 (h 



dX(i 



where in the second inequality we have considered the area inside and outside of [—1, 1] separately. 
Inside it we use the following estimate which follows from Lemma [9] (whenever {x, y} C K, with 
S :=y-x, 0^ \x\ > \6\): 



(12) 



\e lx — ix — e v + iy\ < \5x\ 1 + 



+ 



8x 2 



= \6 X \ + -6 +- - 



o 



and here the first inequality uses \J\ + z < 1 + z/2 in turn (z > 0). We apply (12) to x = pu, 
y = ps in order to estimate the difference inside [—1,1] in (|11|). 



Now, combining ([9]), (10) and 11 the desired claim follows, at least under scheme 1. 
[Under scheme 2, Q remains unaltered. (10) becomes: 



(13) 



\n h gh(p)\ < \v f ^A([-l, l]\[-A/2, h/2}) + [ \u\dX(u)) 

1 \ y[-i,i]\[-$,|] / 



Finally in the estimate (11) we apply (12) to x = up and y = p(s + /isgn(s)) so that we only have 



\x\ > i<5 and hence \e tx — ix — e %v + iy\ < 5\x\ + |<5 2 — (11) now reads: 



V c h 3 {e i{s+hsm(s))p -l)-ipn h - f (e lpu - 1- ipul [_!,!] (w))dA(«) 

» \ rm JR\\-h/2,h/2] 



(14) 



= 6Z fe \{0} 



„ 3ft. ., 3ft 2 , 7 ^3/1 



[-l,l]\[-h/2,h/2] 



/2,h/2] 

dX(u) 



By combining ([£]), ( |13[ ) and (14) the result follows also under scheme 2.] 

Finally, if b = oo (under 5° = V = 1; for all h < 2), denoting £(<5) := x 2 dA(x): 



(15) 



h / cos(Tip) - 1 p : 



f-'o 



V ft 2 



< p 4 ft 2 C(ft/2)/4! 
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(16) 



(- — ) d\(u) - f (e lpu ~ 1 - ipuj d\(u)\ 

(sin(pu) — pu) d\(u)\ 



p 2 u 2 

cos(pit) — 1 H — --r— ) d\(u 
h h-i \ 2! 

' 2 ' 2 I v 

< p 4 (ft/2) 2 £(ft/2)/4! + |p| 3 (ft/2)£(ft/2)/3! 



(17) 



\v h 9h(p)\ = 



h ( sin (ftp) 



\u\d\(u] 



S~] c h 3 {e lsp - 1) - ipfi h - / (e lp " - 1 - ipul[_ M] (u))dA(it) 

» > rm JK\[-fe,4l 



sSZ fc \{0} 



E 



s6Z h \{0} 



- 2piil[_i i i](u) + ipsl[_x 5 i] (zz) rfA(ii) 



(18) < / d\(u)+p 2 ^ f \u\d\(u)+p 2l -(f) A([-1,1]\[-A/2,V2]) + ~«V2) 

where in the second inequality we have now made the following estimate (whenever {x, y} C 
with <5 := y — x, 7^ \x\ > \5\): 

(19) |e la; -e ly -ix + iy\ < \e x - e ly - ix + iy + x 2 /2 - y 2 /2\ + X -\x 2 - y 2 < ~ \S\x 2 + \S\ \x\ + h 2 
using Lemma [9| We apply (fl9|) to x = pu, y = ps. 



[Under scheme 2, (15) and (16) remain unchanged and (17) reads: 



(20) 



M 9h(p)\ < 



ft 2 / 3ft 



■A(t-l,l]\[-ft/2,A/2])+ , 
1 • /r -i.i]\[-#,#] 



\u\d\(u 



Finally, in place of (18) we get: 



cl(e av — 1) — ipfi h — / {e pu — 1 — ipul[—i t i] (u))d\(u) 



sGZ fe \{0} 



(21) < f |P| 



dA(u) +p 2 ^ 

-1,1] 2 J[-l,ll\[-i/2A/2] 



,3ft. 

Y 



where ^4 is some positive constant whose precise value is not important.] Now, combining (15) 



(16), (17) and (18) under scheme 1 (resp. (15), (16), (20) and (21) under scheme 2) yields the 
desired inequality. □ 

The above entails convergence of ^> h (p) to ^(p) as h ! pointwise in p £ R. Levy's continuity 
theorem \Q\ P- 326] and stationarity and independence of increments yield at once Remark pi 



4. Rates of convergence for transition kernels 



Finally let us incorporate the estimates of Proposition 12 into an estimate of D^ T (x,y) (recall 



the notation in ([6])). Assumption [I] and Table [I] are understood as being in effect throughout from 
this point onwards. Recall that \^ h — \&| < (J 2 \fh\ + lAdh\ + \lh\- 
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We assume given a function K with the properties: < K(h) f oo as h J, and if(/i) < \ 
for all sufficiently small h (this condition will always be easily checked and suitable choices 
of K will be identified later). Divide the integration regions in Q and ([5]) into five parts: 
(— oo,— £], (— x, — K{h)), [— K(h), K(h)], (K(h),j^), [^,oo). Separate (via a triangle inequal- 
ity) the integrals in the difference D^ T (x,y) accordingly and use the triangle inequality in the 
second and fourth region, thus (we let s := T — t > 0): 



27vDt T (x,y) < 



+ 



-n/h 



+ 



/h 



tp{*(p)s}|dp + 



-K(h) 



+ 



K(h) 



]exp{<I> h (p)s}\ + |cxp{*(p)s}|) dp + 



K(h) 



cxp{^(p)s} — exp{^I'' l (p)s} dp. 



K(h) 



Now gather the terms with |exp{^I / (p)s}| in the integrand and use \e z 
estimate the integral over [—K(h), K(h)], to arrive at: 



II < - 1 (z e 



to 



(22) 
(23) 
(24) 





' f-K(h) 




< 


J + 






J — oo 


)K(h) 




" f-K(h) 


r- 


+ 


J-t + 








)K(h) 



exp{fy(p)s}\dp + 
ex.p{fy h (p)s}\dp + 



+ 



f K{h) 

I | exp{^(p)s}| I exp{s V (p) - V(p) } - 1 ) dp. 

J-K(h) V ' 



For convenience we let (22), resp. (23) and (24) stand for the first, resp. second and third term on 
the RHS of the above inequality. We deal with each of them separately; the first two will always 
die off quickly, governed by the rapid decay of the exponential function. By a suitable choice of 
K(h) we shall be able to achieve supr^M w^yi — ^| — > as h ]. 0, while still maintaining 
sufficient decay in the other terms. Then for all sufficiently small h, s\^ h — \P| < 1 everywhere on 
[-K(h),K(h)]. Write: 



e «l* h -*l_l = „|¥*_*| + £ 

k=2 



{s\f h - *|)* 
k\ 



< s|*' l -*| + (s|* /l -*|) 2 e s l* h -*l < s|* /l -*| + e(s|*' l -*|) 2 



Manifestly the second term will always decay strictly faster than the first (so long as they are not 
0). We shall use this estimate repeatedly in what follows. 

Suppose first a 2 > 0. We have then | exp{^(p)s}| = exp{ — ^a 2 sp 2 } everywhere in p G R and 



exp{$> h (p)s}\ = exp <^ a 2 s 



cos(hp) 



h 2 



< 



exp 



- 2 

-a sap 



the latter inequality holding for all p G [— ? , ?] and some a > 0. Now, with this estimate as the 
integrand, extend integration in (23) to infinity (i.e. to the whole of K(h), K{h))). 

Similarly, when a 2 = 0, then, as in Remark [ZJ for some C > and all large enough p, 
| exp{*(p)s}| < exp{— Cs|p| e }. Also, the real part of ^f h (p) is majorized on [— |, — K(h)]L)[K(h), |] 
by (for some ft > 0): 
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(25) 



E 



c h s (s + sgn(s)hf 



s€Z h \{0},\s\+h< ^ 



4|p| 



< 



u 2 d\{u) < 



-Pp 2 < — Pp 



< -c'\p\ 



7tt 
4|p| 



3h 7tt_ 
2 >4|p| 



3hi 
2 



u 2 d\(u) 



4 1 p '4 p J 

for some C > C > 0, where we have taken into account that for all small enough h, we will have 
K(h) large enough, so that ^ is small enough for the limit inferior in Assumption jlj to produce 
A > (if 2 A is the limit inferior under question, then there exists some 5 > 0, such that for 
all < r < 6, J^_ rr ^u 2 dX(u) > Ar 2 ~ e ). Hence \ exp{^ h (p)s}\ < exp{-C" s\p\ e } the latter for all 
p £ [— ¥, — lif(/i)] U [-RT(/i), It follows that, with this estimate as the integrand, we can safely 
extend the integral in ( |23| ) to infinity (i.e. to the whole of K(h), K(h))). 

The rest of the analysis relies heavily on the choice of K{h) in each particular instance. In the 
interest of clarity, we shall therefore deal with the individual cases one by one. What follows is 
summarized by Theorem [2] 

Notation-wise we let DCT stand for Lebesgue dominated convergence theorem. 



4.1. Brownian motion with drift: a 2 > 0, A(R) = 0. Let K(h) := y^l^log^. It is then 
manifest that @ and Q will be 0(h 3 ). As for ([24]), notice that K(h) n h m -)■ as h | 0, 
whenever {m,n} C N. Hence for all sufficiently small h, we have s\^> h — ^\ < 1 everywhere on 
[—K(h),K(h)]. Using the inequality e s l*' I -' I 'l - 1 < s\^ h - *| + e(s\^ h - ^|) 2 , the second term 



contributes a term which is 0{h A ) to (24), since exp{ — la 2 sp 2 }dp integrates every polynomial in p 



absolutely (use the estimates of Proposition 12). Finally we have: 



fK(h) poo -i / I 13 „4\ 

J ^ ^ \e W {*(P>}\\*\p) ~ *(p)\sdp <h 2 J exp{--a 2 sp 2 }s [\^ + a 2P -j dp. 



Performing the integration we arrive at (26) of: 



Proposition 13. When a 2 > and c = 0, we have for all sufficiently small h (depending on \i, 
a 2 and s ): 



(26) 



A s (h) < h 2 



1 Id 1 

i + 



1 



3vr fi 4 s 8\/27r {sa 2 ) 3 / 2 



+ 0(h 3 



Hence A s (h) = 0(h 2 ). With a 2 s = 1 and \i = we have limsup^ A s (h)/h 2 > l/(8\/2vr), proving 
that in general the convergence rate is no better than quadratic. 

Proof. It remains for us to establish the last assertion. Take x = y = 0. Notice that if one 
can write D^ T (0,0) > A(h) — B{h) with lim/40 B(h)/h 2 = 0, then it is sufficient to maintain 
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limsup^iQ A(h)/h > 1/(8v2tt). Careful analysis of our estimates shows that we can take: 

1 j-K{h) x 



A(h) = -H / e W {-y}h(p)dp\. 
2n J_ kw 2 

By DCT, we have A(h)/h 2 — > ^ exp{— ^p 2 }p 4: / 4\dp and the claim follows. □ 



Remark 14 (RD). In general, suppose we seek to prove that f > converges to no faster than 
g > 0, i.e. that limsup^ f(h)/g(h) > C > for some C. If one can show f(h) > A(h) —B(h) and 
B = o(g), then to show limsup fe | f(h)/g(h) > C, it is sufficient to establish limsup fe | A(h) / g(h) > 
C . We refer to this extremely useful principle as reduction by domination (henceforth RD). 

4.2. Compound Poisson process with diffusion: a 2 > 0, < A(R) < oo. Let K{h) : = 
yj log \. It is then manifest that (22) and (23) will be 0(h 2 ). Finally for (24), following a 
similar argument to that in the previous case, apart from a term which is 0(h 2 ) we have: 

K{h) , , ,„. , 4 , , has 



K(h) 



hcs 1°° 1 
exp{ty(p)s}\\l h (p)\sdp < — J exp{--a 2 sp 2 }\p\dp. 



Performing the integration we arrive at (27) of: 

Proposition 15. When a 2 > and < c < oo, we have for all sufficiently small h (depending on 
fj,, a 2 , s and c): 



(27) A s (h) <h^^ + 0(h 2 " 



Hence A s (h) = 0(h). With a 2 = s = 1, fi = and X = 5(^1/2+^-1/2) one has lim sup h ^ A s (h)/h > 
showing that convergence in general is indeed no better than linear. 

Proof. With regard to the last assertion, via RD by a careful analysis of our estimates (again 
considering x = y = 0) , one need only observe for the sequence h n = I as n — > 00 the lim sup as 
n — > 00 of j^tH f-KCh ') ex P{ — \p 2 ~ 1+ cos (p/2)}^/v„(p)^p|- By our careful of choice of A and then h n 
it is clear that l hn (p) = 2^(e ip ^ 2 - hn ^ - e ip / 2 ) = cos(p/2)(cos(p/i„/2) - l) + sin(p/2) sm(ph n /2). 
It is therefore sufficient to observe: 



1 /-AUM 1 
limsup- — — I / exp{ — -p — 1 + cos(p/2)} sin(ph n /2) sin(p/2)dp|. 

ra->oo 27T/l n J-K(h n ) 2 

By DCT it is equal to: 

1 Z" 00 1 
I := — / psin(p/2) exp{--p 2 - 1 + cos(p/2)}<ip. 
2vr J Q 2 

The numerical value of this integral is (to one decimal place in units of 2ir) 0.4/(27r), but 
we can show that / > analytically. Indeed the integrand is positive on [0,6]. Hence 
el > sin(l/2)e cos ( 3 / 2 ) jf pe- p ^ 2 dp - e J 6 °° pe~ pi 'I 2 dp = sin(l/2)e cos ( 3 / 2 )[ e - 1 / 2 - e" 9 / 2 ] - e" 17 . Now 
use sin(l/2) > (1/2) • (2/tt) (which follows from the concavity of sin on [0,7r/2]), so that, very 
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crudely: / > (l/7r) e - 1 / 2 (l - e " 4 ) - e~ 17 > (l/^) e - 1 / 2 (l/2) - e~ 17 > (l/e 2 )e-^ 2 (l/e) - e' 17 > 



-17 



> 0. 



□ 



4.3. a 2 > 0, k(0) < oo = A(R) . Let K(h) := W ^log£. (22) and (23) are then <3(/i 2 ). As for 
(24), modulo a term which is o(h) (following the discussion above) we have it bounded by: 



hs 
~2 



exp{— -a 2 sp 2 } (d\p\ + bp 2 ) dp. 



Hence (28) of: 



Proposition 16. When a 2 > 0, and b < oo = c, we have for all sufficiently small h (depending 
on n, a 2 , s, d and b): 

bs 



(28) 



A s (h) < h 



1 d 1 

+ 



2ira 2 1^(o 2 sfl 2 _ 
In general convergence is no better than linear. 



o(h). 



Proof. To show the validity of the final remark, we can for the most part recycle the example from 
the previous subsection. Indeed just take x = y = 0, h n = l/3 n , a 2 = s = 1, /j, = and 



A = - 



1 1 °° 

^0*3/2 + £-3/2) + ^ J^l/3 fe + S ~l 



/3 k 



k=l 



Because the second term in A lives on U ne pj<S/i„ it is seen quickly (via RD) that one need only 

consider (to within non-zero multiplicative constants): 

rK(h n ) 1 1 00 

limsup / — sm(ph n /2) sin(3p/2) exp{--p 2 + (cos(3p/2) - 1) + Y7cos(p/3 fc ) - I)} dp. 

n^oo J-K(h n ) h n 2 fr[ 

By DCT it is sufficient to observe that: 

/•27T/3 1 p 2 00 1 f 00 1 

/ sin(3p/2)pexp{--p 2 + (cos(3p/2) - 1) - — ~^) d V ~ / pexp{--p 2 }cip > 0. 

To see the latter, note that the second integral is immediate and equal to: e -1 / 2 ^ 71 "/ 3 ) 2 . As for the 
first one, make the change of variables u = 3p/2. Thus we need to establish that: 

f'TT 

A := (4/(9e)) / sin(u)uexp{-« 2 /4 + cos(u)}du - e~ 1/2( - 2n/3)2 > 



Next note that — u 2 /4 + cos(w) is decreasing on [0, ir] and the integrand in A is positive. It follows 
that 



A> — 
~ 9e 



J itsin(u) exp | — ~ (J^) + cos ( ^ | du + J usm(u) exp ^ — ~ + cos^^j-du 



-2tt 2 /9 



Using integration by parts, it is now clear that this expression is algebraic over the rationals in 
e, \/3 and the values of the exponential function at rational multiples of ir 2 . Since this explicit 
expression can be estimated from below by a positive quantity, one can check that A > 0. □ 
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4.4. a 2 > 0, k(0) = oo. With regard to the choice of K(h) we have to be a bit more careful here: 



K(h) :-- 



i 



1/4 



C(At/2) y ' Then (22) and (23) will decay strictly faster than Q(h/2) while at the same 
time we still have s\^> h — < 1 everywhere on [—K(h),K(h)] for all sufficiently small h. Thus all 



our other considerations remain the same and we have (22) bounded from above for all sufficiently 
small h and modulo a term which is o(£(/i/2)), by: 



s ( C(h/2) + - 7 (ty2) 



exp{ — -o~ 2 sp 2 }p 2 dp. 



Proposition 17. When a 2 > 0, b 
characteristic triplet and s): 

1 



oo ; we have for all sufficiently small h (depending on the 



(29) 



A s (h) < 



2^(aW~ 2 (c(V2) + ^(V2))+«(C(/»/2l). 



In general the convergence is no better than C,(h/2). 



Proof. To show that convergence is no better than Q(h/2) in general, consider again x = y = 0, 
2 - s = 1, /i = 0, /i„ = l/3 n and A = X^i^fcOW + 8- Xk ) where x„ = \h n and w„ = l/x n 



a 



(n G N). Notice that oi := /[_i,i]\[_/ ltl /2,h„/2] « 2 ^(u)=2 ££ =1 and ^ := Y, s &s hn \{o} c s n s 2 
2 ^2k=i( x k~ ^2 L ) 2 ' w k an d A := (Ti — (72 = 2£(/t n /2)— 7 (/t n /2) > Q{h n /2). It is sufficient to show then, 
that /Jjfe \P 2 exp{^(p)}dp converges to a strictly positive value as n — > oo, which is transparent 
(since ^ is real valued). 

4.5. a 2 = 0, «(0) < oo = A(R). Taking 2T(/t) := 



□ 



2-2 

C" 



t log r we see that ( 22 ) and ( 23 ) will be 0{h 



Finally we have (24). Note that on account of our choice of K, we will yet again have s\^ h — < 1 



everywhere on [-K(h), K(h)] for all small enough h. Hence, as usual, we arrive at the first assertion 
of: 



Proposition 18. When a 2 
no better than linear. 



and b < oo = c we have A. s (h) = 0(h). In general convergence is 



Proof. To establish convergence is no better than linear in general, take s = 1, x = y = 0, /i = 
and A = | X^feLi w k{$x k + $-x k ), where x n = h n = l/2 n and w n = 2 n / 2 (n G N). By checking 
Orey's condition on the decreasing sequence h n , Assumption [T] is satisfied with e = 1/2 and we 
have b < oo = c. fi h = by symmetry. 

Going through our estimates and using 



RD, 



we 



need only show that 



limsup^ jjjj |5j[m exp{ty(p)}(ty h (p) — ^(p))dp\ > 0. Moreover, using the estimates of 
Proposition |12| and again RD, it suffices to show: 



lim sup — 

n— >oo tin 



/K(h„) I /• 

exp{*(p)} Yl / 



-K(hn) 



(cos((s + sgn(s)h n )p) — cos(pu)) d\(u) dp 



seS hn \{oy 



> 0. 
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and via RD get rid of the first two terms (i.e. they contribute to B rather than A in Remark 14). 
It follows that it is sufficient to observe: 



lim sup — 



r K(h n ) 00 / n \ 

/ exp{ V(cos(p/2 fc ) - l)2 fe / 2 } V 2 k ' 2 sin(p/2 fc ) h n pdp 
J-K(K) k=1 Vfc=i / 



> 0. 



Finally, via DCT and evenness of the integrand, we need only have: 

POO / 00 \ OO 

/ ^2 fc / 2 sin(p/2 fc ) pexp{]T(cos(p/2 fc ) - l)2 fc / 2 }dp ^ 0. 
Jo \k=i J k=i 

One can in fact check that the integrand is positive, as the following lemma shows, and thus the 
proof is complete. □ 

Lemma 19. Let for I E N U {0} and p > 0, ip l (p) := Y^T=i 2fc/2 s'm(p/2 k ). Then the ip l are all 
positive on (0, 00). 

Proof. Since ?/>' +1 (2p) = ^/2tfj l (p), it is sufficient to show the claim for 1/j =: ip. The key observation 
for ip is the identity ip(2p) = sin(2p) + y/2i/j(p), which we use as follows. 
Clearly (i) ^1(0,1] > an d (ii) for p E (1, 2] we have: 

V>(p) = 4>{2?-) = sin(p) + V2tp(p/2) > sin(l) + \/2>(l/2) = V(l) > 5/2 > 1/(^2 - 1) =: A , 



since ip is strictly increasing on [0, 1], sin2 > sinl, 5/2 > l/(v2 — 1) and > 5/2 (see below). 
We now claim that if for some B > we have tp\(p,B] > an d ^{(b^b] > A)> then V'l(o,2_B] > 
and ■0I(2_B,4B] > ^0) and hence the assertion of the lemma will follow at once (by applying the 
latter first to B = 1, then B = 2 and so on). So let 2p E (25,45], i.e. p E (B,2B\. Then 
^>(2p) = sin(2p) + v^G?) > -1 + V2/{V2 - 1) = 1/(^2 - 1), as required. 

As a technical aside we should verify the slightly non-trivial observation ip{V) > 5/2. Indeed, 
?/>(l) is, by definition, a series of positive terms. Hence, by using sin(l/2 fc ) > sin(l)/2 fc for k > 0, 
ip(l) > sin(l) + sin(l) ^^ =1 2~ fc / 2 = sin(l)[l + ^ - ]. Next, the series expansion of sin around 
and the properties of alternating series yield sin(l) > 1 — 1/6 + 1/5! — 1/5! = 5/6. Thus 
^(1) > |[1 + "-/=-[] > 5/2, as is easily verified. □ 



4.6. a 1 = 0, k(0) = 00. Take K(h) := ^(^/ 2 ) ) • Then (22) and (23) decay strictly faster than 
((h/2). As for (24), we still have s\^ h — ^\ < 1 everywhere on [— K(h), K(h)] for all small enough 
h. Hence the first assertion of: 
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Proposition 20. When a 2 = 0, b = oo, we have A s (h) = 0(C(h/2)). In general, convergence is 
no better than £(/t/2). 

Proof. To establish the second assertion, take s = 1, x = y = 0, fi = 0, hk = Xk = l/2 fc , Wk = l/%k 
(k G N) and A = Ylh=l w k(o~x k + S- Xk ). Checking Orey's condition on the decreasing sequence Xk, 
Assumption [l] obtains at once with e = 1. Observe that [i h = and \E' is real valued. 
By a careful analysis of our estimates, using RD, it is sufficient to observe that: 



lim sup 

n— too 



K(h n ) p 2 exp r^(p)} 

-K(h„) 



X{h n /2) 



E 



s + /i n sgn(s)) 2 



s&S hn \{0} 



u 2 d\(u) 



dp 



> 0. 



Direct calculation yields the expression in the square brackets as being equal to: 

n n n n 

^(Xk + K) 2 W k - ^ x l w k = 2/t n E XkWk + h n E Wk = 2 (^(^n/ 2 ) + 7(W 2 )) ' 
k=l k=l k=l k=l 

Hence it is sufficient to observe that J^Z^p 2 exp{^(p)}dp > and this follows, since ^> is real- 
valued. □ 
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